Exploring Symmetry Breaking at the Dicke Quantum Phase Transition 



K. Baumann, R. Mottl, F. BrenneckeQ and T. Esslinger 
Institute for Quantum Electronics, ETH Zurich, 8093 Zurich, Switzerland 

We study symmetry breaking at the Dicke quantum phase transition by coupUng a motional 
degree of freedom of a Bose-Einstein condensate to the field of an optical cavity. Using an optical 
heterodyne detection scheme we observe symmetry breaking in real-time and distinguish the two 
superradiant phases. We explore the process of symmetry breaking in the presence of a small 
symmetry-breaking field, and study its dependence on the rate at which the critical point is crossed. 
Coherent switching between the two ordered phases is demonstrated. 
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Spontaneous symmetry breaking at a phase transition 
is a fundamental concept in physics [1 . At zero temper- 
ature, it is caused by the appearance of two or more de- 
generate ground states in the Hamiltonian. As a result 
of fluctuations, a macroscopic system evolves into one 
particular ground state which does not possess the same 
symmetry as the Hamiltonian. Finding a clean testing 
ground to experimentally study the process of symmetry 
breaking is notoriously difficult as external fluctuations 
and asymmetries have to be minimized or controlled. The 
protected environment of atomic quantum gas experi- 
ments and the increasing control over these systems offer 
new prospects to experimentally approach the concept 
of symmetry breaking. Recently, rapid quenches across a 
phase transition were studied in multi-component Bose- 
Einstein condensates [T-T and optical lattices E]. 
Such a non-adiabatic quench causes a response of the 
system at correspondingly high energies. Therefore, a 
central characteristic of a phase transition, which is its di- 
verging susceptibility to perturbations, remains partially 
hidden. 

In this work we study the symmetry breaking process 
while slowly varying a control parameter several times 
across a zero-temperature phase transition. Compared 
to quenching, this allows us to explore the low energy 
spectrum of the system which probes its symmetry most 
sensitively. For very slow crossing speeds we identify the 
presence of a residual symmetry breaking field of vary- 
ing strength. Larger values of this residual field can be 
correlated to the repeated observation of one particularly 
ordered state. For increasingly steeper ramps across the 
phase transition the influence of the symmetry breaking 
field almost vanishes. 

We investigate the symmetry breaking in the motional 
degree of freedom of a Bose-Einstein condensate coupled 
to a single mode of an optical cavity. Our system real- 
izes the Dicke model frWS. which exhibits a second-order 
zero-temperature phase transition p^OHIS] . The broken 
symmetry is associated with the formation of one of two 
identical atomic density waves, which are shifted by half 
an optical wavelength [HI El HH [15] . Using an interfero- 
metric heterodyne technique, we monitor the symmetry- 



breaking process in real time while crossing the transi- 
tion point. A similar technique has been used to test 
self-organization in a classical ensemble of laser-cooled 
atoms [15 , where the symmetric phase is stabilized by 
thermal energy rather than kinetic energy [16 . 

The Dicke model [7 considers the interaction between 
N two-level atoms and the quantized field of a single- 
mode cavity, which is described by the Hamiltonian 
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H = huJoJz + hujCL^a H — i={cL^ + ci)Jx. 
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Here, a and denote the annihilation and creation op- 
erators for the cavity mode at frequency cj, and J = 
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FIG. 1. (Color online) (a) Experimental setup. A Bose- 
Einstein condensate is placed inside an optical cavity and 
driven by a far-detuned standing- wave laser field (wavelength 
Ap) along the z-axis. Phase and amplitude of the intracav- 
ity field are measured with a balanced heterodyne setup (PD: 
photodiodes). (b) Steady-state order parameter (J^) as a 
function of coupling strength A, with corresponding atomic 
density distributions (1.-3.). The order parameter vanishes 
in the normal phase (1.) and bifurcates at the critical point 
Acr, where a discrete Ap/2-spatial symmetry is broken. The 
two emergent superradiant phases (2. and 3.) can be distin- 
guished via the relative time-phase 0. 



{Jx^Jy^Jz) describes the atomic ensemble with transi- 
tion frequency cjq in terms of a pseudospin of length N/2. 
The cavity light field couples with coupling strength A to 
the collective atomic dipole Jx- In the thermodynamic 
limit, the Dicke model exhibits a zero-temperature phase 
transition from a normal to a superradiant phase when 
the control parameter A exceeds a critical value given by 
Acr = sJojujQ 1 2 ^1"QUT2] . Simultaneously, the parity sym- 
metry of the Dicke Hamiltonian, given by the invariance 
under the transformation (a^Jx) (— a, — J^^), is spon- 
taneously broken [13 . While parity is conserved in the 
normal phase with (a) = = (Jx)^ two equivalent su- 
perradiant phases (denoted by even and odd) emerge for 
A > Acr, which are characterized by ( J^^) ^ and (a) ^0, 
respectively (Fig. [T]3). 

In our experiment [8 we couple motional degrees of 
freedom of a Bose-Einstein condensate (BEC) with a 
single cavity mode using a transverse coupling laser 
(Fig. [l^). Within a two-mode momentum expansion of 
the matter- wave field, the Hamiltonian dynamics of this 
system is described by the Dicke model (Eq. [l]) [3 [H [17] 
where the effective atomic transition frequency is given 
by ujQ = 2uJr with the recoil frequency Ur = hk'^/2m, the 
atomic mass m and the wavelength Ap = 27r//c of the 
coupling laser. The frequency and power of this laser 
controls the effective mode frequency uj and the coupling 
strength A, respectively [8]. Above a critical laser power, 
the discrete Ap/2-spatial symmetry, defined by the opti- 
cal mode structure u{x,z) = cos{kx) cos {kz), is sponta- 
neously broken and the condensate exhibits either of two 
density waves (Fig.[T]3). Correspondingly, the atomic or- 
der parameter (J^^), given by the population difference 
between the even {u{x, z) > 0) and odd {u{x^ z) < 0) 
sublattice, exhibits a negative or positive macroscopic 
value, while the emergent coherent cavity field oscillates 
(for UJ ^ hi) either in {(j) = 0) or out of phase {<p = tt) 
with the coupling laser. 

As described previously [8 , we prepare BECs of typi- 
cally 2 X 10^ ^^Rb atoms in a crossed-beam dipole trap 
centered inside an ultrahigh-finesse optical Fabry-Perot 
cavity, which has a length of 176 jam. The transverse cou- 
pling laser at wavelength Ap = 784.5 nm is red-detuned 
by typically ten cavity linewidths 2k. = 27r x 2.5 MHz from 
a TEMoo cavity mode, realizing the dispersive regime 
UJ ^ ujQ of the Dicke model. We monitor amplitude and 
phase of the intracavity field in real-time using a bal- 
anced heterodyne detection scheme (Fig. [l^). Due to 
slow residual drifts of the differential path length of our 
heterodyne setup, which translate into drifts of the de- 
tected phase signal of about O.Itt/s, we cannot relate 
the phase signals between consecutive experimental runs 
separated by 60 s. 

To observe symmetry breaking, we gradually in- 
crease the coupling laser power across the critical point 
(Fig.[2|L). The transition from the normal to the superra- 
diant phase is marked by a sharp increase of the mean in- 
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FIG. 2. (Color online) Observation of symmetry breaking and 
steady- state superradiance. Shown are simultaneous traces 
of (a) the coupling laser power P, (b) the mean intracavity 
photon number riph, and (c) the relative time-phase between 
coupling laser and cavity field (both averaged over 150 /is). 
The coupling laser frequency is red-detuned by 31.3(2) MHz 
from the empty cavity resonance and the atom number is 
2.3(5) X 10^. Residual atom loss causes a slight decrease of 
the cavity photon number in the superradiant phase. 



tracavity photon number (Fig. |2]3). Simultaneously, the 
time-phase between the two light fields locks to a con- 
stant value, implying that the symmetry of the system 
has been broken (Fig. The observation of a con- 

stant time-phase above threshold confirms that the sys- 
tem reaches a steady-state superradiant phase in which 
the induced cavity field oscillates at the coupling laser 
frequency. When lowering the laser power to zero again, 
the system recovers its initial symmetry and a pure BEC 
is retrieved, as was inferred from absorption imaging af- 
ter free ballistic expansion. 

To identify the two different superradiant states 
(Fig. [T]3), we cross the phase transition multiple times 
within one experimental run (Fig. [3^). Above threshold, 
the corresponding phase signal takes always one of two 
constant values. From multiple traces of this type we 
extract a time-phase difference of 1.00(2) x tt between 
the two superradiant phases, where the statistical error 
can be attributed to residual phase drifts of our detection 
system. 

If the system was perfectly symmetric, the two ordered 
phases would be realized with equal probabilities, when 
repeatedly crossing the phase transition. However, the 
presence of any symmetry-breaking field will always drive 
the system into the same particularly ordered state when 
adiabatically crossing the critical point. We experimen- 
tally quantify the even-odd imbalance by performing 156 
experimental runs (similar to Fig. |3^) , in each of which 
the system enters the superradiant phase ten times within 
Is. A measure for the even-odd imbalance is given by 
the parameter e = (mi — 777,2) /lO, where 777,2 ^ ^1 de- 



note the number of occurrences of the two superradiant 
configurations in individual traces. In 73 % of the traces, 
the system realized ten times the same time-phase, corre- 
sponding to the maximum imbalance of e = 1 (Fig. [sJd). 
However, 12 % of the runs exhibited an imbalance below 
0.5, which is not compatible with a constant even-odd 
asymmetry. 

We attribute our observations to the finite spatial 
extension of the atomic cloud. This can result, even 
for zero coupling A, in a small, but finite population 
difference between the even and odd sublattice, deter- 
mined by the spatial overlap O between the atomic col- 
umn density n{x^z) (normalized to N) and the opti- 
cal mode profile u{x,z). This asymmetry enters the 
two-mode description (Eq. [T]) via the symmetry-breaking 
term 2h\0{a) + a)/\/N ^ and renormalizes the order pa- 
rameter {Jx) by the additive constant O. The resulting 
coherent cavity field below threshold drives the system 
dominantly into either of the two superradiant phases, 
depending on the sign oi O. In the experiment, the re- 
sulting even-odd imbalance is likely to change between 
experimental runs, as the overlap integral O depends Ap- 
periodically on the relative position between the mode 
structure u{x^ z) and the center of the trapped atomic 
cloud, with amplitude Oq. We can exclude a drift of the 
relative trap position by more than half a wavelength Ap 
on the timescale given by our probing time of 1 s, as it 
would lead to equal probabilities of the two phases, pre- 
tending spontaneous symmetry breaking. 

The openness of the system provides us with direct 
experimental access to the symmetry-breaking field pro- 
portional to O. Indeed, we detect a small coherent cavity 
field (nph < 0.02) in the normal phase whose magnitude 
varies between experimental runs. In all runs exhibiting 
an imbalance of e = 1 (Fig. [sj)), the relative time-phases 
of the cavity field detected below and above threshold are 
equal. Furthermore, the even-odd imbalance increases 
significantly with the light level observed below thresh- 
old. Post-selection of those 10 % of the runs with the 
smallest light level yields a much smaller mean imbal- 
ance (Fig. [sJd, inset). 

In general, the influence of a symmetry breaking field 
becomes negligible, if the mean value of the order param- 
eter, induced by this field, is smaller than the quantum or 
thermal fiuctuations present in the system. From a mean- 
field calculation performed in the Thomas- Fermi limit for 
A/' = 2 X 10^ harmonically trapped atoms, we estimate 
a maximum order parameter of Oq = 40 for zero cou- 
pling strength, corresponding to an even-odd population 
difference of 40 atoms. This value is much smaller than 
the uncertainty AJ^^ = \fN /2 = 224, given by vacuum 
fiuctuations of the excited momentum mode. Therefore, 
one expects in the extreme case of a sudden quench of 
the coupling strength beyond Acr, that the apparent sym- 
metry is spontaneously broken, resulting in nearly equal 
probabilities of the two superradiant phases. 
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FIG. 3. (Color online) (a) Cavity time-phase (red, averaged 
over 30 jis) for a single run, and corresponding time sequence 
of the coupling laser power P (dashed), (b) Probability dis- 
tribution of the imbalance e (see text) for 156 runs, where the 
phase transition was crossed at a rate of A/ Acr = 18(3) s~"^. 
The inset displays the distribution of post-selected data (see 
text), (c) Mean imbalance (dots) as a function of the rate 
A/ Acr at which the transition was crossed (extracted from 
356 runs in total), and theoretical model (solid line). The 
error bars indicate the standard error of the mean of e and 
systematic changes of Acr during probing. 



In the experiment we determined the even-odd imbal- 
ance e for increasingly larger rates A/ Acr at which the 
critical point was crossed, i.e. in an increasingly non- 
adiabatic situation (Fig.jsj^). As the transition is crossed 
faster, the mean imbalance between the two superradiant 
phases decreases significantly and approaches the value 
e ~ 0.25 corresponding to the balanced situation (Fig.jsj^, 
dashed line) . This indicates that the effect of the symme- 
try breaking term can be overcome by non-adiabatically 
crossing the phase transition. 

Our observations (Fig. ^) are in quantitative agree- 
ment with a simple model based on the adiabaticity con- 
dition known from the Kibble-Zurek theory [TSllQ . We 
divide the evolution of the system during the increase of 
the transverse laser power into a quasi-adiabatic regime, 
where the system follows the change of the control pa- 
rameter, and an impulse regime, where the system is ef- 
fectively frozen. After crossing the critical point, fiuctu- 
ations of the order parameter, which are present at the 
instance of freezing, become instable and are amplified. 
The coupling strength which separates the two regimes 
is determined by Zurek's equation [18 |C/CI — with 
C = (-^cr — A) /Acr and the energy gap between ground 
and first excited state given by A = hwQ^l-X^/Xl, for 
a; >cJo [ElIIIl. 

We deduce the probability with which the system 
chooses the even phase, Peven = /o^P(©)dO, from the 
probability distribution p(0) at the instance of freezing, 
where O denotes the shifted dipole operator O = Jx-\-0. 
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FIG. 4. (Color online) Coherent switching between the two 
ordered phases. After adiabatically preparing the system in 
one of the two superradiant phases, the coupling field is turned 
off for a time r. Displayed is the steady- state cavity time- 
phase A0 (averaged over 0.5 ms) after turning on the coupling 
field, referenced to the value recorded before the turning-off. 
Each data point corresponds to a single measurement. The 
dashed line shows the time evolution as expected from the 
two-mode model. 



In the thermodynamic limit the distribution p{Q) be- 
comes Gaussian with a mean value (6) = (Jx) + O 
and a width determined by the quantum fluctuations of 
the order parameter AJx. These values are determined 
from the linear quantum Langevin equations based on 
the Dicke model [17 including the symmetry breaking 
term. Besides the decay of the cavity field we also take 
into account dissipation of the excited momentum state 
at a rate 7 = 27r x 0.6 kHz. This value was deduced from 
independent measurements of the cavity output field be- 
low threshold [20] . 

From the steady-state solution of the quantum 
Langevin equations we find that the mean order parame- 
ter (©) grows faster in A than its fluctuations. If O > 12 
the order parameter exceeds its uncertainty already be- 
low critical coupling. Thermal fluctuations are neglected 
in this analysis. For our typical condensate temperatures 
of about 100 nK quantum fluctuations dominate as long 
as C > 0.005. We account for shot-to-shot fluctuations of 
the overlap O by suitably averaging over the position of 
the harmonic trap. The solid line in Fig.jsj) shows a least 
square fit of our model to the data with the single free 
parameter Oq. We obtain a value of Oq = 77 which is 
in reasonable agreement with the theoretically expected 
value of Oq =40. This verifies the predominance of the 
considered symmetry breaking field over other possible 
noise terms. 

Finally, we experimentally demonstrate coherent 
switching between the two ordered states. To this end we 
suddenly turn off the coupling laser field after adiabati- 
cally preparing the system in one of the two superradiant 
phases. The atoms are then allowed to freely evolve ac- 
cording to their momentum state occupation, giving rise 
to standing-wave oscillations of the atomic density dis- 
tribution. In the two-mode description this corresponds 
to harmonic oscillations of the order parameter (J^) at 
frequency 2co'r. We probe this time evolution by turning 



on the coupling laser after a variable oflP-time r, thereby 
deterministically re-trapping the atoms either in the ini- 
tial or in the opposite superradiant state. As expected, 
we observe regular 7r-jumps in the difference between 
the steady-state phase signals measured before and af- 
ter the free evolution, with a frequency of 2cJr (Fig. [4j 
dashed line). The inertia of the atoms traveling at finite 
momentum causes the 7r-jumps in Fig. 4 to occur before 
those times at which the order parameter has evolved by 
an odd number of quarter periods. 

In conclusion, we have experimentally monitored sym- 
metry breaking in the Dicke quantum phase transition 
and identified the interplay between a residual symme- 
try breaking field, fluctuations and the crossing speed. 
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